








X 



ANDREEV REFLECTION 
AND THE SEMICLASSICAL BOGOLIUBOV-DE GENNES HAMILTONIAN: 

RESONANT STATES 



A. BENSOUISSI i'2, N. M'HADBI ^ & M. ROULEUX 



^ Universite du Sud Toulon- Var, and Centre de Physique Theorique 

CPT, Case 907, 13288 Marseille Cedex 9, France 

^ Universite de Tunis El-Manar, Departement de Mathematiques, 1091 Tunis, Tunisia 

o. 

^S| . Abstract: We present a semi-classical analysis of the opening of superchannels in gated mesoscopic 



SNS junctions. For perfect junctions (i.e. hard- wall potential), this was considered by [ChLeBl] in 
the framework of scattering matrices. Here we allow for imperfections in the junction, so that the 
complex order parameter continues as a smooth function, which is a constant in the superconducting 
banks, and vanishes rapidly inside the lead. We obtain quantization rules for resonant Andreev states 



Q-i' near energy E close to the Fermi level, including the determination of the resonance width. 



0. Introduction. 

Bogoliubov-de Gennes Hamiltonian is a 2 x 2 matrix V{x,^) defined for (x,^) € T*R, which 
describes the dynamics of a pair of quasi-particles (hole/electron) in a 1-D metallic lead connecting 
^.^ ^ 2 superconducting contacts. Diagonal terms are of the form =b(^^ — fi{x)), where /x(x) stands for 

cn I the chemical potential, while the off-diagonal interaction with the supraconducting bulk is modeled 

through a complex potential, or superconducting gap, Aoe*"^*'^ at the boundary ; due to the finite 



^^ I range of the junction, we may consider that the interaction continues to a smooth function x i— > 

L^ ! A{x)e^'^^^''^ on a neighborhood of the lead (say, the interval [—L,L]. ) So we assume that A(x), 

we will call henceforth the "gap function", is a smooth positive function, increasing on x > 0, with 
A(x) = 0, |x| < xi < L and A(x) = Aq, |x| > X2 > L (ignoring the fact that A shows typically isolated 
zeroes (vortices) in the supraconducting bank). In the same way, we will assume that (/)(x) = sgn{x)(j) 
takes only 2 values. The chemical potential /i(x) will be extended also to a smooth positive function 
on a neighborhood of [—L,L], constant and > Aq for |x| > X2- As is usual for a metal, we assume 
that // and A are even in x, which provides this model with the CPT symmetry. 

We introduce a "Planck constant" h > 0, which stands for the ratio of L to the characteristic de 
Broglie wave-length, and take usual h-Weyl quantization, 

(0 I) Vix hD)-( ^^^-^' - ^(^) A(x)e^*W/2 

^^■^> ^ ^^' ^^^^ - V A(x)e-*<^(-)/2 -{hD^f + /i(x) 

An electron e~ moving in the metallic lead with energy < E' < A (measured with respect to Fermi 
level Ep) and kinetic energy Kj^{x) = /i(x) + yjE"^ — A(x)^ is reflected back from the supraconductor 
as a hole e"*", with kinetic energy K-{x) = fi{x) — \J E"^ ~ A(x)^, injecting a Cooper pair into the 
bulk. When inf [.^^^j /i(x) > E, and 7^ 0, this process yields so called phase-sensitive Andreev states, 
carrying supercurrents proportional to the (/)-derivative of the eigen-energies Ek{h) of 'P(x, hD^)- Since 
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V{x,hDx) is self-adjoint, there is of course also an electron moving to the left, and a hole moving 
to the right (in fact, V{x,hDx) is the Hamiltonian for 2 pairs of quasi-particles), for no net transfer 
of charge can occur through the lead in absence of thermalisation. So we stress that Bogoliubov- 
de Gennes Hamiltonian is only a simplified model for superconductivity, and that a more thorough 
treatment should also take into account the self-consistency relations coupling the quasi-particle with 
the gap function A(x) and the phase 0(x), that we treat here as "effective potentials" (see [KeSo]). 

In the case where A(x) is a "hard-wall" potential, this was studied in [ChLesBl], [CaMo] in 
the framework of scattering matrices. In [BelfaRo] , we derived semi-classical quantization rules for 
Andreev states near energy E^ from a microlocal study of the Hamiltonian in the "inner region" 
A(x) < E alone. For simplicity, we assumed that A(x) varies linearly near E, namely if xq G]a;i,a;2[ 
is such that A(xo) = E, then /i(x) = fi = Const, and A(x) = E + a{x — xq) near xq. 

Here we want to take also into account the "outer region" A(x) > E (i.e. |x| > xq) of the 
junction, entering the supraconducting bulk. As a matter of fact, the microlocal solutions, purely 
oscillating in A(x) < E, acquire a complex phase in A(x) > E, which is of course related to phase- 
space tunneling. We make the assumption that the junction is extended, in such a way that the 
quasi-particle turns into a resonant state before creating a new Cooper pair, its dynamics still being 
governed by Bogoliubov-de Gennes Hamiltonian. So we assume that /i(x) and A(x) are defined on the 
entire real line, taking constant values for |x| > X2 > L, so that (0.1) can be defined as a self-adjoint 
operator on L^(R) (8) C^. We will translate the usual theory of analytic dilations [ReSi] in the context 
of CPT symmetry, and find semi-classical resonances near a "scattering" Andreev level, i.e. complex 
correction to the real eigen-energies Ek{h) of V{x, hD^)- 

1) The real part of the resonances. 

The bicharacteristic set in {^ > 0} at energy E, of the form detV{x,S,) — E = 0, or ^^ = K±{x), 
consists of : (1) two real curves p± over [— xq, Xq], joining smoothly to a close curve at the "branching 
points" a' = {—xo,S,o) and a = (xo,^o) (so to make p^ U p- diffeomorphic to S-*^) ; (2) complex 
branches p^ over ] — oo,— xq], and p-^ over [xo,+oo[ respectively. They all have a vertical tangent 
at a, a'. We complete this picture by reflection on the x axis, denoting the corresponding branching 
points by b' , b. 

a) Microlocal solutions supported on p± 

First we recall from [3] the construction of distributions microlocalized on the Lagrangians p±, 
and verifying the FT symmetries of the problem. We denote the parity operator by ^ : u{x) — t- u{—x), 
and the time reversal operator by I : u{x) — t- u{x). 

Definition 1.1: We call "admissible C'^-valued Lagrangian distribution" an oscillatory integral 
(1.1) /(5,v?)(x,/i) = (27r/i)-^/2 /■ e*'^(^'®''^)/'*5(x,e;/i)de 

with the following properties : (1) {p{x,Q,h) denotes a non degenerate phase-function, and 

S{x, 6; h) = So{x, 0; h) + hSi{x, 6; /i) + • • • 



a C^ -valued amplitude (i.e. a classical symbol in h), Sq = {^ y ) possibly depending on h (with the 
property that (/»(x) = sgn{x)(j)) ; (2) The symbols X = X{x, Q,h),Y = Y{x, 0, h) have their principal 
part \yj ~ -^(^' ®i ^)(y'); A € C, proportional to a real vector (y?), depending also on {x,Q;h). 

Of course, all these functions may depend on additional parameters. One of the main problem 
consists in finding microlocal solutions near the branching points a, a'. Due to PT symmetry, it 
suffices to focus on o = (xo,^o)- In /i-Fourier representation, the Hamiltonian takes the form 

^^ - fM e^'t'/^iE-ahD^-axo) 

'{E — ahD^ — axo) — ^^ + /f^ 



V^{^hD^,0-\ ^-.^/2. 



where ^u = ^q i^ ^ constant, equal to the value of the chemical potential at xq. Consider the equation 
(V°'{—hD^,^) — E)U = 0, where C/ = I ^ I . Clearly, the system decouples, and to account for time- 
reversal symmetry, it is convenient to introduce the scaling parameter /3 = ^/a{2^Q)~^''^ , together with 
the changes of variables ^ = ^o(±2/3^' + l). The functions M±/3(f) = (^2_^_^)-i/2g-i(£;-axo)C/afe^2 
satisfy a second order ODE of the form 
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(1-2) {P±f^{-hD^',i',h)-^)u^p{i' 

with El = {2(o)~'^E, and 



P±^{-hD^,,e,h) = {hD^,f + (e' ± /3e 



/2 



,2 



^^'^^ + h\2^^)-^p\2p^e ± 2/3C' + \ + Er) [P^^ + PC - i^i)"' 

Operators Pfj and P-fj are unitarily equivalent, and so have the same spectrum. Up to the O^h"^) 
term, P±fj{—hD^i,^',h) have the structure of an "anharmonic oscillator", with "potential wells" at 
^' = 0,^1/(3 separated by a "barrier" at ^' = =Fl/(2/S). It is also well known [HeSj] that, viewed 
as a /i-PDO of order 0, microlocally defined near (x',^') = 0, P = P±i3 can be taken to the normal 
form of a harmonic oscillator, away from the "barrier". More precisely, there exists a real-valued 

CxD 

analytic symbol F{t,h) = F±^{t,h) ~ ^Fj{t)h^, defined for t € neigh(O), Fo(0) = 0, F(,{0) = i, 

3=0 

Fi{t) = Const., and (formally) unitary FIO's A = A±fj whose canonical transformations ka defined 
in a neighborhood of (0,0), are close to identity and map this point onto itself, such that 

A*F{P, h)A = Po = ^{{hD^f + v^-h) 

Define the large parameter v by F{ ^^ ^ , h) = vh. So u = u-t solves (1.2) microlocally near (0,0) 
iff t; = A*u solves Weber equation (Pq — i'h)v = microlocally near (0,0), when uh ^ -^-^ is small 
enough. The well known parabolic cylinder functions D^^ and D_jy_i, provide with a basis of solutions 
of h{{hDy)'^ + rf — h)v = vv. We shall use £)_,y_i, and write 



:i.4) V = A*u^p = Y^ a^-^-'^D.,_i{ie{h/2)-^'^n) 



(-U-1) 

•Jt 
e=±l 



for complex constants a^- ^a • 

These microlocal solutions can be expressed in the spatial representation by taking inverse h- 
Fourier transformation ; expanding integrals of the type (1.4) by stationary phase, both pieces of 
bicharacteristics p-i- contribute to Ui.^±j3 near a. Microlocal solutions near a' are deduced by PT 
symmetry. 

Once microlocal solutions U'^'I'^~ have been obtained that way near the branching point a, it 
is standard to extend them up to a' as WKB solutions {U^'p'^~ )ext, taking advantage that V has 
simple characteristics away from a, a'. When A(x) = 0, i.e. for —xi < x < xi, they are completely 
decoupled, which means that the solution is either a pure electronic state, i.e. colinear to the vector 
(„) of C^, or pure hole state, i.e. colinear to (-^) . Otherwise, they are a superposition of electronic/hole 
states. We summarize these constructions in the : 

Proposition 1.2: For x < xq near xq, there are 2 basis of oscillating microlocal solutions of {V"' — 
E)U = indexed by e = ±1 : 



u;':,±,i^,h'), Y.K:^:±~A^^^') 



p = ±l /9=±1 



Here the branch with p = p± = ±1 is microlocalized on p±, i.e. the part on p_|_ (^ > ^o near a), 
belongs to the electron state, while the part p- (^ < ^q near a) belongs to the hole state ; they satisfy, 
for p = ± : 

(1.5) u;--y = u;'-y + o{h) 

and 

(1-6) ui'-y = [/»;;;_-; 

Each of these solutions is an admissible C^ -valued lagrangian distribution in the sense of Definition 
1.1. Divide all microlocal solutions by the trivial factor e«'^/4g«-Eo5o//i ^ j^^ — j^ _ q^xq. Then with the 
notations of (1-4) the general solution of (V^ — E)U = is of the form 

p,£ 

The solutions near a' are given by symmetry, e.g. U"^ ^~_^^ = ^XC/"'^ T . Moreover both microlocal 
families can be extended as WKB solutions (satisfying Definition 1.1) along the bicharacteristics. 

Note also that in this region where /i(x) is a constant, Up^e,±fi = e^^^"^^l^p,e,±fi;h' with Up^e,±ii;h' 
oscillating on a frequency scale 1/h' = l/{ah), so if we think of the slope a to be large, Up^e,±i3 
behaves as a plane wave e^^^"'^ , modulated by a slow varying function. 

b) Real holonomy and approximate Bohr-Sommerfeld quantization condition. 

The microlocal kernel Kh{E) oiV — E on] — xo,xo[xR_(_ can be viewed as a 4-D fibre vector 
bundle J-h{E) of admissible Lagrangian distributions over S^. We characterize the real part of the 



resonances as the set of Ek{h) near E such that this fibre bundle is trivial. We start by computing 
the holonomy of that bundle. 

First we normalize the basis in K}^{E) obtained in Proposition 1.2 using generalized Wronskians 
introduced in [HeSj], [Ro]. Namely, let x = x" be a smooth cut-off supported on a sufficiently small 
neighborhood of a, equal to 1 near a, uj± = wj a small neighborhood of p± n supp[P,x"]) and 
Xcj± = Xu± a cut-off equal to 1 near uj±. We take Weyl /i-quantization of these symbols, and for 
C/,FeKft(^), wecall 

y\^.^{u,v) = {x^J-[VMU\V) = {x<.^l[P,x]u\v) 

the microlocal Wronskian of {U,V) in u}±. This is a sesquilinear form on Kh{E), and Wuj±{U,U) is 
independent, modulo error terms 0{h°°), of the choices of x^ and xZ± as above. Taking into account 
both contributions of p± we define also 

W{U, V) = W^+ {U, V) + W^_ {U, V) 

For each microlocal solution U = U'^'^a~ , it turns out that W{U,U) have asymptotic expansions in 

h', of the form Wo (£;,/3)-F/i'u>i (£;,/?) H , with wo(£;,/3) > 0. 

Given x = X"; l^t now x = x** be a new cut-off equal to 1 on the support of x"; and to outside 
a slightly larger set. For U,V £ Kh{E) we set {U\V)~ = {xU\V). Then it is easy to see that there is 
an orthonormal basis of Kh{E) for the "scalar product" {U\V)~, which is at the same time orthogonal 
for yV{U, V) (everything being defined modulo 0{h°°) .) This allows to find Ve = V^'^^~ of the form 
(1.7) such that (K|K,)~ = ,5e,e', {e,e' = ±l),W{V±,V±) > 0, and W(y+,y_) = 0. Of course, by the 
symmetry '^IV{x,hDx) = V{x,hDj;)X^ , such normalized microlocal solutions exist as well near a'. 
The Lagrangian distributions 

eld =Xc.--f^[P,x]u,;^ 

and similarly F^\^'^~ span the microlocal co-kernel K^{E) oiV — Ein] — xo,xo[xR_|_, as e = ±1. 
The same holds for or G'^'7'^~ obtained by replacing U'^'a''~ by the "orthonormal basis" V^'r^~ 
as above. 

Because of Proposition 1.2, the normalized microlocal solutions V^ ^J" are related to the ex- 
tension of the normalized microlocal solutions V^'r^~ along the bicharacteristics by a monodromy 

matrix M"'" = I -^"^ }^^ ) G t/(2). Similarly, we obtain M" '" by extending from the left to the 

V"21 "22/ 

right, and due to symmetry, M°- '°- = {M°-'°- )~^ = {M°-'°- )*. Diagonal entries of these matrices are 
given by action integrals along p± (see e.g. [Ro]). Off-diagonal terms are 0{h') and can be computed 
with the help of the Wronskian (in the ordinary sense) associated with the system (V — E)U = (see 
e.g. [Ba]). 

The quantization condition is satisfied, precisely when the rank of that system drops of one 
unit (actually, because of degeneracy, of 2 units), i.e. when dim Kh{E) = dim. Kl{E) = 2. This 
amounts to set to zero the determinant of some Gram matrix Gram(i?, h) expressed in the basis 



(^e"«~'^~ ■•^e' '-H~ ) ^"^^ i^t'p"' '^£''-fl~ )• So i? = Ek{h) is ail eigenvalue, modulo 0{h°°), of 
V{x, hDx), corresponding to an Andreev state, iff det Grain(£', h) = 0. 

Here we note the sensitivity of the energy levels Ek{h) with respect to (j). In the "hard- wall" limit 
a — 7- oo, we recover the quasi-particle spectrum, of the form cos^ = cos( ^ ^ — 2arccos( ^ )) 
for some smooth function g (see [CayMon], [ChLesBl]). 

2) The imaginary part of the resonances. 

The considerations above are not sufficient to account for exponentially small corrections to 
Ek{h). Further information will be extracted from a Grusin problem. 

a) Microlocal solutions with complex phase. 

Microlocal solutions, computed in the real phase space, are purely oscillating in the metallic 
lead [— Xo,a;o]. To get information in the "superconducting part of the junction", we need use "in- 
finitesimal" invariance by time reversal and conjugation of charge. The substitutions ^' i— > zizi^' , or 
equivalently, /? i— ?• =bi/3, leave invariant equation (1.1), with a new operator P±i;3 in Fourier-Laplace 
representation. Microlocal solutions of {P±ii3 — E)u = are constructed similarly and, on the real 
domain, independently of those of {P±ii — E)u = 0. 

Thus the fibre bundle of microlocal solutions on R x R_|_ (i.e. microlocal kernel oiV — E) splits 
as F^{E)®Fh{E)®F>{E), where we recall Fh{E) from Sect.l, and F^'^{E) are 2-D (trivial) fibre 
bundles over R. 

Nevertheless, taking advantage that the coefficients are analytic near a, a', there is a way to 
couple Ff^ ' (E) with Fh{E) in the complex domain. This, together with the assignment that the 
global section be "outgoing" at infinity, accounts for complex holonomy. 

First we investigate complex holonomy near a, and consider the family of operators, obtained 
by extending P±ji along a path {e*'*'/3, < 7 < 27r} in the complex plane ; similarly, we consider the 
family of Lagrangian distributions obtained by extending u^(^') along that path. They will solve 
(1.1) iff 7 = 0, ±7r/2, vr. 

These solutions are related through their Lagrangian manifolds as follows : consider (for sim- 
plicity) the principal part of P-t^ and P±i^., namely Qfj{—hD^',^') = {hD^iY + {^' + /3^'^) and 
Qip{-hD^,,i') = {hD^.f + (C' + i/3^'^)^ The potentials being equal for ^' = and C' = -2/(1 + i)/3, 
the real Lagrangian manifold />+ near a extends analytically along the loop {e*'''/? : 7 G [0, 27r]} in the 
complex domain, so that it intersects p^ at —2/(1 -|- i)j3 for 7 = 7r/2. We can argue similarly for the 
other branches. Actually, both p± and p^ are branches of a single 2-sheeted Riemann surface, with 
complex "turning points" . 

We can assign to this analytic manifold microlocal solutions for Peiip as in (1.2) with complex 
phase, which yields in turn solutions oi {V — E)U = Q for relevant values 0, ib^,7r of the parameter 
7 ; these solutions are very similar to the C^"-|-«'s given in Proposition 1.2. The monodromy operator, 
acting on microlocal solutions, is known as connection isomorphism^ see [DeDiPh] and references 
therein, and also [Fe], or [Ro,Sect.4,g] in the case of a system. This connection isomorphism is given 
by a matrix N"' G U{2), whose entries are expressed in term of exponentials of action integrals 
computed along Stokes lines between the complex turning points. 



Let us consider next the conditions at infinity : for |x| > X2, V has constant coefficients, so 
we make an analytic dilation of the form x i— t- exp[(sgnx)^] x, i) > 0. Plane waves with positive 
momentum have the phase exp[ix(^i it 1^2) /h] where ^1 it i^2 = (a*o =t i\/^Q — i?^) , according to 
the choice of p±''' ■ Analytic distorsion is turned on for |x| large enough, and 1? in the complex upper- 
half plane. We denote by V^ the distorted operator. So for all Im t? > small enough, we can make 
the "electronic state" (resp. "hole state") exponentially decaying at +00 (resp. —00), which models 
the scattering process e"*" — t- e~, and similarly for the scattering process e~ — ?■ e"*", thus preserving 
conservation of charge. 

b) A Grusin problem and the w^idth of resonances. 

Following a classical procedure in Fredholm theory, we can translate the original eigenvalue prob- 
lem for V into a finite dimensional problem via the Grusin operator [HeSj3,Sect4] ; this is essentially 
the isomorphism (^^(R) (g) C^)/KhiE) -> Ran(P - E) C L^(R) (g) C^. Here KhiE) denotes the 6-D 
microlocal kernel of V — E inH x R_(_, restricted to the set of outgoing functions defined above. For 
V = V^, we consider G{E) = 0(1!}, E) of the form : 



R+ 

(2.1) ' 6 

R.{xi,---,xe) = Y.x,G„ i?+[/=(([/|G,))^<^.<g 

where the Gj's range over the basis of co-kernel K^E) consisting of G" ^, G" _^, G" ifj,G°_ _^^ (or 
their analytic continuation at the branching points). 

At this point we make the following remark : Since our Grusin operator (2.1) involves only 
positive frequencies, it cannot be associated with the self-adjoint operator Vi) (for real d-). But 
resonances are due precisely to a breaking of time-reversal symmetry, and their imaginary part is 
computed by introducing a /i-Pseudo-differential cutoff $(x,/iD^) supported in {^ > 0}. Because 
negative frequencies will be eventually removed, we may best think of (2.1) as a short-hand notation 
for the "full" Grusin operator Q{E), that would take into account the negative frequencies as well. 

For all /i > small enough, Q{E) is bijective, with bounded inverse 

/ £o{E) £+{E)' 
^^^'-\ £_{E) £-+{E) 

and has the property, that E is an eigenvalue of V iff det£ \-{E) = 0. The construction of £{E) is 

carried as in [HeSj], [Ro], selecting solutions according to the prescriptions above. Matrix £ \-{E) 

decouples modulo 0{h°°), with a 4 x 4 block conjugated to Gram(i?) ; the interaction with the 
"incoming hole" and "outgoing electron" occurs through the "turning points" in the complex domain, 
involving the connection isomorphisms N°',N"' . 

For complex (i?, E), we note that (V^ — E^ = V-^—E. Applying distorsion to the Grusin operator 
as well, we get : 



We can prove that Gi'd, E) is well-posed for all i? G C small enough, with inverse £{^, E). Recall from 
[Ro, Prop. 7.1] the following identity : 

Proposition 2.1: Let <& G C^(R^;R). With the notations above 

[R-{^,EY^£+{^,E)V£-+{'&,E) -{[R-{^,Er^£+{^,E)X £.+{{}, E))* 
= £+{^,EY[V^M£+{^,E) 

In the self-adjoint case, the corresponding statement would be ^\R*_£jY)*£ h is self-adjoint". The 

determination of the width of resonances then goes as in [Ro], though it is somewhat more com- 
plicated due to the structure of £-+{"&, E). Take W{i},E) € Kex£_j^{'d,E), and set A{i},E) = 
[R-{i},E)*^£+(J),E)\* . From (2.2) and the identity 

{W{'d,E)\A{'d,E)£_+{§,E)W{'d,E)) ~ {A{^,E)£_+{^,E)W{'d,E)\W{'d,E)) =0 

we get 

{A{^,E)£_+{d,E)W{{},E)\W{'d,E)) - {A{^,E)£_+{^,E)W{'d,E)\W{^,E)) 



(2.3) 



{£+ii&,Ey[P^,<^]£+{i&,E)Wii^,E)\Wii},E)) 



Evaluating both members of this equality gives an implicit equation for the imaginary part of the 
resonance, showing that behaves like exp[— 2 J ^ dx/h'], where r C C is a path connecting the complex 
branching points in p^ f] p.^.. 
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